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We consider a dynamo wave in the solar convective shell for the kinematic aw-dynamo 
model. The spectrum and eigenfunctions of the corresponding equations are derived an- 
alytically with the aid of the WKB method. Our main aim here is to investigate the 
dynamo wave behavior in the vicinity of the solar poles. Explicit expressions for the inci- 
dent and reflected waves are obtained. The reflected wave is shown to be relatively weak 
in comparison to the incident wave. The phase shifts and the ratio of amplitudes of the 
two waves are found. 



1 Introduction 

The solar cycle is a well-known manifestation of the magnetic activity of 
the Sun. The physical mechanism responsible for this activity is thought 
to be the dynamo generation of the large-scale solar magnetic field (see e.g. 
Zel'dovich, Ruzmaikin and Sokoloff, 1983). There are a number of numerical 
models of the solar and stellar dynamo waves which reproduce many features 
of the observed solar magnetic activity (Brandenburg, 1994; Riidiger and 
Brandenburg, 1995; Tobias, 1997 etc.). 

It is important to support these numerical investigations by approximate 
analytical solutions of the dynamo equations. Recently, new analytical meth- 
ods of solving the mean field dynamo equations have been developed for the 
case of a very intensive generation. However, analytical solutions were ob- 
tained only in the simplest cases. Kuzanyan and Sokoloff (1995) derived 
the asymptotic solution of the Parker migratory dynamo equations (Parker, 
1955). Meunier et al. (1997) and Bassom et al. (1997) have considered 
asymptotic solutions of the Parker migratory dynamo in the nonlinear regime. 

These results, as well as the Parker's equations themselves, have a limited 
domain of applicability. In particular, they do not reflect the details of the 
dynamo wave behavior in the very vicinity of the poles. In the present paper, 
we consider dynamo effects in the framework of a more general model, taking 
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into account the convective shell curvature. In the limiting case of very large 
dynamo numbers, we derive consistently the spectrum and eigenfunctions of 
the dynamo equations and describe the dynamo wave near the solar pole. 

According to the modern observational results, the dynamo wave and the 
wave of sunspots, observed in the butterfly-diagrams, propagate equator- 
wards within the main spatial domain (i.e. far from the solar equator and 
solar poles). Besides, there is a weak dynamo wave in the vicinity of the 
solar poles which was observed by Makarov and Sivaraman (1983). This 
wave propagates pole-wards. In this paper, we explore this polar dynamo 
wave by solving the generalized equations of the Parker migratory dynamo. 
We show that the incident dynamo wave reflects from the pole. The reflected 
wave is relatively weak compared to the incident wave. We demonstrate also 
that the process of dynamo wave reflection reduces to a phase jump and the 
reflection is not accompanied by any sharp changes of the wave amplitude. 

Our paper is structured as follows: In Sec. 2, we derive the equations in 
the Parker's model from the general equations of the mean field electrody- 
namics. In Sec. 3, we solve these equations asymptotically with the aid of the 
WKB method, using the fact that the dynamo number is very large for the 
Sun. We obtain the spectrum and the eigenfunctions describing the dynamo 
wave far from the pole. In Sec. 4, we reduce the equations to a simplified 
form valid in the vicinity of a pole. In Sec. 5, we solve these asymptotic 
equations and show that our solution describes two waves: the incident wave 
and the reflected wave. In Sec. 6, we match two asymptotics found in Sees. 
3 and 5 and derive the ratio of the amplitudes of the incident and reflected 
waves and the corresponding phase shifts. In summary section, we discuss 
our results and make some simple numerical estimates. 

2 Basic Equations 

Here we obtain a generalization of the well-known Parker migratory dynamo 
equations which takes into account the curvature of the convective shall. The 
large-scale magnetic field generation in a turbulent flow of a differentially 
rotating electrically conducting fluid is governed by the following equation 
(see Krause and Radler, 1980): 
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— = Vx («B) + Vx (v x B) + /3AB, 



(1) 



where B and v are the large-scale (mean) magnetic and velocity fields cor- 
respondingly, a is the helicity coefficient and (3 is the turbulent diffusivity. 

Looking for a kinematic axisymmetric eigensolution of Eq.(0), we present 
the magnetic field as a superposition of the poloidal and toroidal fields as 
follows: 

B(r,t) = [B p (r)+B f (r)]e^, (2) 

where B t = ^0, 0, B\ is the azimuthal component of the magnetic field, 
Bp = i?rot (0, 0, A) is the radial component of the magnetic field (R is the 



solar radius) and 7 is the eigenvalue to be found. Thus, Re 7 is the magnetic 
field growth rate and 27r(Im7)~ 1 is the dynamo wave period. 

In the new terms, Eq.(|l|) takes on the following dimensionless form: 
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where we have introduced the following dimensionless constants: 
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Here, G 



v ' " * s ^ ne ra dial gradient of the mean angular velocity and 

G' = y^q~ * s ^ s meridional gradient. Note, that lengthes are measured in 
units of the solar radius, time in units of the diffusion time r^f = R 2 /(3, a 
and G in units of their maximal values. 

We suppose that the differential rotation is more intensive than the mean 
helicity and that G(r, 9) weakly depends on the latitude. Thus, the terms 
in Eq.(§) which contain R a and G' can be omitted. This approximation is 
known as aw-dynamo model. 
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We also suppose that the process of intensive magnetic field genera- 
tion takes place in a thin shell . Parker (1955) has noted, that in this 
case, two-dimensional dynamo equations could be reduced to effectively one- 
dimensional equations by their averaging over r (see also Kuzanyan and 
Sokoloff, 1996). Below, we explore the corresponding one-dimensional prob- 
lem. 

All the assumptions, made above, yield a significant simplification of 
Eqs.([| f|) and we have: 



1 A = a{6)B + j- e 
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Here D = R a Ruj is the dimensionless dynamo number, which characterizes 
the intensity of the generation sources, function A(6) = R~ 1 (A(9, r)) is pro- 
portional to the averaged azimuthal component of the vector-potential and 
B{9) = (B(9,r)) is the averaged azimuthal component of the magnetic field 
( (...) means averaging over r), 9 is the latitude measured from the solar 
pole. Helicity coefficient a{6) is also averaged over the shell section. We 
suppose that a(0) ^ 0. Note, that after averaging, the diffusive terms take 
the following form (see also Proctor & Spiegel, 1991): 
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and they can be taken into account by a redefinition of eigenvalue 7. 

Mention, that a typical estimate of the dynamo-number for the Sun is 
\D\ ~ — 10 4 or even more. Below, we investigate Eqs.(|5|, ||) analytically 
taking into account the large value of the dynamo number. We derive the 
asymptotical spectrum corresponding to these equations and the asymptotic 
behavior of the dynamo wave. We consider two spatial domains: The first 
one is the main domain which is arranged far both from the poles and from 
the equator and the second one is the domain near a solar pole (9 <C 1). 

Let us note that the simplest approximate form of the one-dimensional 
dynamo equations, was obtained phenomenologically by Parker in 1955 (see 
also Stix, 1989). Those equations, known as Parker's equations follow from 
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equation system (|5], |6[). They appear as a first approximation for the case of 
very short dynamo waves: 



"< A = a[6)B + IF' (7) 

7B = DGW sin^ + ^. (8) 

Let us mention that using the Parker's equations one can obtain the correct 
expressions for eigenvalues 7 in the leading approximation. However, the 
eigenfunctions do not coincide, even in the main approximation, with the 
corresponding solution of Eqs.(|5]||). Note also that equations similar to (|7|, 
H) were solved numerically in the nonlinear regime by Jennings (1991) in the 
main domain. 



3 Asymptotic solution in the main domain 

To obtain the asymptotic solution of Eqs.(|5], [(J) in the main domain, we use 
the WKB approach (see e.g. Maslov and Fedorjuk, 1981). Similar method 
was applied by Kuzanyan and Sokoloff (1995) to solve the Parker's equations. 
In this section, we follow their treatment to obtain the solution of the more 
general equations under discussion. 

Let us rewrite Eqs.(||, |6]) in the spectral form explicitely: 
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where H is a linear differential operator: 
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Let us present the eigenvector and eigenvalues in the following form: 
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where e — \D\ (1/3) , 
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7 = ^ro + -ri + -... (12) 
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Here fi, fix, v, v x and S are complex functions of the latitude and T , Fx 
are complex constants. Emphasize that value e = \D\~a is the true param- 
eter of our asymptotic expansion and its smallness is the main condition of 
applicability of our solution. 

To find the spectrum of Eq.(^), we have to formulate boundary conditions. 
Worledge et al. (1997) have emphasized the role of these conditions. They 
showed that the solution in the linear regime is very sensitive to the changes 
in the boundary conditions (see Tobias et al., 1997 for the nonlinear regime). 
Our asymptotic WKB expansion is not applicable near the boundary, i.e. 
in the very vicinity of the poles and the solar equator. Thus, we can not 
formulate boundary conditions in the explicit form in the framework of the 
asymptotic theory. However, we can match the asymptotic WKB solution 
applicable in the main domain with the solution applicable near the pole. 
The boundary condition for the latter is that the magnetic field is limited 
everywhere including the pole. To perform the matching one must require 
that the asymptotic solution decays for 9 — > and 9 — > |. Otherwise, the 
field grows exponentially in the vicinity of the boundary and the matching 
is impossible. The other requirement is that the asymptotic solution is a 
smooth function. These two conditions are found to be sufficient to obtain 
an asymptotic spectrum. 

Note, that the dynamo number can have any sign. The sign determines 
the direction of the dynamo wave propagation. We choose the most realistic 
case in our solution. Using the observational fact that the dynamo wave 
propagates equator-wards in the main domain, we accept the negative sign 
for D that corresponds to this behavior. 

Equations of the first approximation can be obtained by substituting 
Eqs.(|TT]) and (0) into Eq.(||) and equating the terms of the minimal power 
of e. This yields: 



r + k{9f -a{9) 
iG(9)sm9k(9) T + k(9) 2 




0. (13) 



A non-trivial solution of this equation exists only if the determinant of the 
matrix in the left-hand side of Eq.([U|) vanishes. This condition leads to 
so-called Hamilton- Jakobi equation (we use the standard terminology of the 
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Figure 1 : Roots of Hamilton- Jakobi equation (|l4|) in the complex plane k plotted as 
functions of latitude 9 for the chosen value of IV The points of matching of different 
branches are circled and the corresponding values of 9 are given. 

WKB-method, see e.g. Landau and Lifshitz, 1958): 

(r + k 2 f + iak = 0, (14) 

where k{6) = is the generalized momentum (by analogy with the semi- 
classical approximation in quantum mechanics). We have also introduced a 
new function a(6) = G(6)a(8) sin8. 

Note that the Hamilton- Jakobi equation in our problem coincides with 
the one for the Parker's equations. Its solution is described in the paper 
of Kuzanyan and Sokoloff (1995). Here, we will recall shortly the reason- 
ings that allow to calculate spectral parameter To from the Hamilton- Jakobi 
equation. 

The accepted decay condition (the solution should be small near the 
boundary) can be rewritten in the terms of the generalized momenta as 
follows: 

Im/cl^^r > 0, lmk\ e ^ < 0. 
These expressions play the role of boundary conditions. 
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Equation ( |HD is the algebraic equation of the fourth order and thus it 
possesses four branches of roots. However, none of these branches satisfies 
the boundary conditions. It means that we should construct our solution by 
matching two (or more) branches. To explore this case, we denote 



H(k,9) = (T + k 2 f + ika(9). 

In terms of function H(k, 9), it is easy to formulate the conditions of crossing 
of two branches at a point 9': 

H[k(9'),9'] = 0; 

(15) 

^[k(9'),9'} = 0. 
These equations can be solved explicitly and we have 

(1 , 2) 3a(9')h^ n(3) _ 3a(9') 2 s 

1 — 78 5 1 — « • 
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(3) 

As one can see, the real part of T is negative, which corresponds to a 
decaying solution. The two other conjugated values correspond to growing 
solutions of our interest. The leading mode is the most important from the 
physical point of view. Thus, we choose 9' to be the point where function 
a(9) reaches its maximum (see Fig.l). It was shown (Galitski and Sokoloff, 
1998), that only this choice of T leads to smooth eigenfunctions of the Parker 
equations. 



From Eq. (|T3|) , we have: 




r + k^ 

— ik sin 9 , 



a(9), (16) 



where a (9) is to be defined from the equation of the second approximation. 
To obtain this equation, we substitute Eqs. flTTD and (|l|) into Eq.@ and 
equate the terms of the first power of e. This yields: 

( T + k 2 -a [ik' + ikctg9 -T l } f i + 2ikfi' \ 

\ -iksin9 T + k 2 ) l^J ~ \[ik' + ikctg9 - v + 2ikv' - [fj l sia9)'J ' 
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Here, for the sake of simplicity we neglect the term containing G'(8), sup- 
posing that the angular rotation weakly depends on the latitude in the main 
domain. Note, that ( |T7| ) differs from the corresponding expression for the 
Parker's equations. 

We see that the matrixes in the left-hand sides of Eqs.(|l3|) and ( |P7D are 
identical. As we have seen the corresponding matrix is degenerate. Thus, 
Eq. ( |17| ) possesses solutions only if the Fredholm resolvability condition is 
satisfied, i. e. if the vector in the right-hand side of Eq. ( |17D is ortoganal to the 
eigenvector of the adjoint equation. It yields so-called transport equation: 



r\ - ik' 1 + 
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Using the fact that the function in the brackets in the right-hand side of 
equation ( |18|) vanishes at 0' (recall that 0' is the point at which a is maximal) 
and following the standard procedure of the WKB theory (see Maslov and 
Fedorjuk, 1981) we can find the spectrum. First, let us denote the following 
functions: 

rv 
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and 




P{0) ctg0. 
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We suppose that P(0), Q(9), <j{0) are analytical functions in the vicinity of 
0' . Furthermore, we present all these functions in the form of Taylor series: 



P($) = [P + P 1 (0)(0 - 0')} {0 - 0'), where P = 3ik' 

3ik'(0') 

Q(0) = Qo + Qi(0) (0 - 0') , where Q = 
For a n we have 

a n (0) = (0-0T[C n + C n+1 (0-0')], 



(21) 
(22) 

(23) 
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where n is an integer parameter, which classifies the eigenvalues. Substituting 
(|!D,(!D and dH) into (PD and setting 9 = 9', we have 

(24) 



P n = I\ n - Q . 
From here we can easily obtain the spectrum: 



The explicit expression for o is 



n 



a(9) 



1 



r 1 - ik' i + 



sin 9 



exp 



0,1,2, 



2k 2 



2ik - 



a 



d9 



2(r + k 



(25) 



(26) 



We emphasize that Eqs.(^) and (|26|), in contrast to the Hamilton- Jakobi 
equation, differ from the corresponding expressions obtained by Kuzanyan 
and Sokoloff (1995). Surprisingly, our solution has a simpler structure despite 
the fact that the initial equations in our problem are more complicated than 
the Parker equations. 



4 Equations Near The Pole 

We now proceed to investigate equations ([5], ^) in the second domain near 
the solar pole. In the following, we have two small values: 9 and e. Since 
9 -C 1, we can rewrite equation Eq.(|5D in the following form: 

. . d 2 A 1 dA 1 . , 

7j4 = q . ( 0) B + _ + -_--A (27) 

Here we keep only the terms of the minimal power of 9. 

Taking into account results of Sec. 3, it is reasonable to assume that 
\A\ ~ £ 2 \B\ and each differentiation of the field multiplies it by e^ 1 (formally 
we can write j| ~ jj). Our assumptions will be confirmed by the results 
obtained below. This yields the following estimates of the terms in (|6|): 



B" ~ - 4 \Al 7 S ~ \\A\, ctg9B' ~ ^\A\, -J—B ~ -^\A\, 
e 4 9e 6 smfc* 2 9 2 e 2 
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6 1 
DG sin6 A' ~ —\A\, DG cos OA ~ — U|. 

The terms corresponding to the last two estimates can be neglected and 
Eq.(|6|) reads: 

d 2 B ldB 1 
^ = W + 6 M - ¥ B ' (28) 
Equations ([27|, |2S| ) can be rewritten in the following compact form: 

Jx{x)A(x) = ^-B(x); (29) 

J 1 (x)B(x) = 0, (30) 

where we have introduced a new variable x = yf—jO and a differential oper- 
ator 

j / \ d 2 Id / 1 
dx z x dx V x / 
Let us note that in Eqs.(^, |3(]) 7 plays the role of an external parameter. 
Its value was obtained in Sec. 3. 



5 Dynamo Wave in The Vicinity of The Solar 
Pole. 

In this section, we solve equations Eqs.(^. ^(]) obtained above. Equation (j3C|) 
is the homogeneous Bessel's equation of the first order for B and Eq.(^) is 
the inhomogeneous Bessel's equation of the first order for A. The magnetic 
field must be finite for all x. It means that we should take Bessel function of 
the first order Ji( X^j clS db solution of Eq.(j30|). Hence, 



B(x) = 2C B J l {x) = C B H^\x) + H\ 



r( 2 )( 



(31) 



Here C B is a constant and H[ 



(1,2), 



are Hankel functions of the first type 



and the second type of the first order. Note, that Hankel functions diverge 
at the pole (9 = 0), but their sum is finite. As we shall see below, 



(2) 

describes the incident wave and H{ (x) describes the reflected wave. 
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The general solution of Eq. (p9|) is a sum of general solution of the cor- 
responding homogeneous equation and a particular solution of the inhomo- 
geneous equation. To find particular solution Ai(x), note that the Wrohski 
determinant of Hankel functions is 



H[ 1] (x)H[ 
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Summarizing, we have 
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Two complex constants and in (|34T) are not independent. Really, 
fields A(x) and B(x) represent two different components of a true eigenvec- 
tor for Eq. (|| |^) . The norm of this eigenvector is an arbitrary quantity, but 
its orientation is prescribed. Our asymptotic solution ([33]) must satisfy this 
condition too. To find a connection between Ca and Cb, one should con- 
sider a spatial domain where both asymptotics ( |TT1) and fl34j) are valid. This 
domain is characterized by the following condition: 



e < e < 1 



(35) 



Since |x| 3> 1, we can use the well-known asymptotic form of Hankel 
functions for large arguments: 



H[ 1] (x 




32L 

4 , 



JL fi -<*-¥) 



(36) 
(37) 



12 



i.e. 



B(x,t)^C B \ — 
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Let us emphasize, that constraint (^) is not a condition of applicability of 
Eg. fl34|) , but a condition of applicability of asymptotics ([3£ 



and ( pEop (see 

Sec. 6). These approximate presentations allow us to perform the matching 
in the explicit form. However, the condition (^) can be satisfied only if e is 
extremely small. Otherwise (if it is not possible to find a domain where the 
value 6 is both much lesser than unity and much greater than e), one should 
perform the matching using the general form of our asymptotic solution. 
Recalling that x = y/—^f6, one can define the phase surface as 



,t) 



ReV-To Imr 
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The phase velocity is 
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(39) 



If ^phase < than the wave propagates pole-wards, if u p hase > than the wave 
propagates equator-wards. 

Evaluating v p hase explicitly, we readily conclude that h[ 1 \x) corresponds 
to the incident wave and '{x) corresponds to the reflected wave. 

The asymptotic form of A(x,t) is 



A(x,t) 
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6 Matching of Asymptotics 

One can see from ([39]) that branches k% and (see Fig.l) correspond to the 
incident wave. Now, we have to find a branch of roots ki(6) that corresponds 
to the reflected wave. We expect that the reflected wave, if exists, should 
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be very weak. It means that this wave decays going equator-wards. Thus, 
we must choose branch k2 to represent the reflected wave. Only this branch 
is arranged in the upper part of the complex plane k and so it describes a 
decaying wave. 

Now, we should match a linear combination of asymptotics (|ll|) for the 
incident and reflected waves and asymptotics fl40| , p8|) . For this aim, it is 
necessary to evaluate approximate expressions for k±(0), cr±(9), A±{6) and 
B±{9) for the case 6 1 C 1 (here and further +'s correspond to the incident 
wave and -'s correspond to the reflected wave). From equations ([□])> (Q), 
( |26|) and (^5|) after a rather long but straightforward algebra we have follow- 
ing expressions (for 9 — > 0): 

k + {9) =k - Ke'^Ve, (41) 
k_(9) = -ko + ne^Vd, (42) 

where K and kn 



K = ili~ (a ,,i > fc o= \' e (43) 
a ± {9) = ie*^, (44) 
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Using these expressions, one can obtain the approximate form of the solution 
of (^, |G|) in the main domain: 
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From the other side we have (see (|4T| |38|) ) : 



Xo\/ |Cb 



3tt ' 



where 
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Our aim now is to express constants C A , Cb and C_ through C + . To 
find the corresponding expressions, one should match the amplitudes of the 
incident and reflected waves in Eqs.(f46|) and (|47| ) at a point 9\. After some 
algebra, we have 

C_ = iC+, (48) 
i 

C B = ^%C + , (49) 

The last step is to find the point Q\ at which the matching is being made. 
In solutions (f46[ ) and fl47f ), the phase-dependences are the same. However, 
phase-shifts appear in the higher order approximations (see Eqs.(|44]) and 
( f45"|) for the explicit expressions). We choose Q\ to be the point at which the 
phase-shifts for different asymptotics of the incident wave are the same. It 
yields 

e 1 = 3£lmf+ (si) 

2/tcos^ v ; 

Mention, that our results are in agreement with all the assumptions made 
in Sec. 4. From Eq . (|50|) , we see that A ~ e 2 B, from Eq . (flTD , we obtain 

^ ~ Y and ~ as we have supposed. However, we see that Q x ~ e 



and this result is on the border of applicability of asymptotics ([36], |37|) . 
Strictly speaking, we should use the exact form of Hankel functions to per- 
form the matching, instead of using their asymptotic representations. It 
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could have some influence on the values ([4^ — pT|) , but not on ([28]). For the 
case of simplicity, we avoid here the complification, connected with the exact 
determination of 6\. However, the physical results of our analysis remain ap- 
plicable for any values of 6\. Moreover, we emphasize, that asymptotics ([36], 
|37|) reflect the behavior of Hankel function with a very good accuracy even 
for | a; | ~ 1 (see Fig. 3) and hence, expressions ([49|-|5TD are well applicable as 
well. 



7 Summary And Discussion 

This paper has derived the analytical expressions for the large-scale magnetic 
field in the vicinity of the solar poles. These expressions can be used to 
compare the relative magnitudes of the incident and reflected waves. 

Taking into account the fact that generalized momenta k for the incident 
and reflected waves differ only by sign at the pole [see Eq. (|43|) 1 , we can 
estimate the ratio of the absolute values of the amplitudes of the reflected 
and incident waves at point 9 as follows (for the exact dependence see Fig. 3): 

R{9) « exp (~ \Imk \e\ 

where for k see fl43f) . Obviously, this ratio depends on 9 and, as one can see, 
the farther from the pole, the weaker the reflected wave is compared to the 
incident wave. Let us assume, that the dynamo- number is D = — 10 4 (or, 
equivalently, e ~ 0.05) and the helicity coefficient is a(9) = cos9. Evaluating 
9\ explicitly we have 9\ ~ 2.3°. This means that for 9 < 2.3° the solution 
( |TTD is applicable, for 9 > 2.3° one should use expressions fl3"4|). Let 9 = 10°, 
than we obtain 

R(10°) « 0.03. 

It follows that the reflected wave is about 30 times weaker than the incident 
wave. However, the incident wave itself is rather weak near the pole against 
the background of the dynamo wave in the main domain. For example, the 
incident wave at 9 = 10° is approximately 10 times weaker than the dynamo 
wave in the domain of generation (see Kuzanyan and Sokoloff, 1995). It 
follows, that the reflected wave is 300 times weaker than the main wave. 

Summarizing, we see that the solar magnetic field in the accepted model 
is described with the aid of the three dynamo waves. The most intensive 



16 



me) 



O.J 



o.r 



o.j 



o.J 



o.r 



o.J 



o.r 



o.J 



o.J 



1.0 








20° 



25° 



6 



Figure 2: The ratio of magnitudes of the incident and reflected waves as a function 
of latitude 9. The solid line represents dependence R{0), obtained with use of (|34|). The 
dashed line corresponds to the approximate dependence R(0), calculated with the aid of 
asymptotics ( |36| , [37|). Both curves are extrapolated for large 8. Asterisks represnt some 
values R{6) calculated using the WKB solution ([ll]). The value D = — 10 4 for all curves 
is accepted. 

wave propagates from high latitudes equator-wards. This wave, of course, 
is well-known both from observations and from numerical simulations. The 
second wave propagates in high latitudes pole-wards. It has essentially lesser 
intensity than the first one. However, it has been detected in the observations 
of the solar subpolar flares (see Makarov and Sivaraman, 1983). In this paper 
we predict the existence of another dynamo wave, which reflects from the 
pole and propagates equator-wards. We show that the reflected wave is very 
weak and it decays exponentially propagating from the pole. It is possible 
however that the reflected wave could be detected in the future specialized 
observations of the solar magnetic activity in the subpolar domain. 

Using the solution obtained, we can also estimate the generation threshold 
as a dynamo number, for which the generation of magnetic field commences: 
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Hence, for the leading mode we have (n=0): 



n(cr) 



Rer A 3 2"/» 



Re Tq / air, 



For a(#) = cos(6 l ) we have Z?o ~ 90.5. Note, that this value is about 
two times larger than the critical dynamo number for the Parker's equations. 
Evaluating similar expressions D£ r ' for n = 1, 2, 3, . . . and taking into con- 
sideration existing estimates of D for the Sun, one can conclude that only 
one additional mode (with n — 1) is excited. One should mention, that even 
provided that = 100, we see, that the true parameter of our asymp- 

totic expansion is not very small and Eq.(|52"D should be considered as a lower 
estimation for the critical dynamo number. 

As it was shown in Sec. 3, the value of T depends on the point 9' at 
which the matching of branches k 3 and k± is made. We have chosen this 
point to be the point of maximum of a{9). The general case for the Parker's 
equations was considered by Galitski and Sokoloff (1998). It was shown that 
the corresponding eigensolutions can not be chosen to be smooth functions 
for any other values of 8'. The same result can be obtained for Eqs.(|| |6|). 



It means, that operator H [see (|10|)] does not possess any other discrete 
eigenvalues except the ones obtained above (p5|). 

It is interesting also to compare our asymptotic theory with so-called 
Maximally- Efficient Generation Approach (MEGA), see Ruzmaikin et al. 
(1990). The latter suggests that the generated magnetic field is localized in 
a small vicinity of the point at which the magnetic field generation sources 
have maximum. Far from this point, the magnetic field appears due to the 
diffusion mechanism only. As one can see, the point of maximum of our so- 
lution is shifted from the point of maximum of function a(9) (this property 
of the solution was emphasized by Kuzanyan and Sokoloff, 1995). It means, 
that in our case the MEGA does not work quite well and one should be 
careful in application of this approach. 

Here, we considered dynamo equations in the kinematic approximation. 
Let us note that in the recent paper of Meunier et al. (1997), it was shown 
that in some nonlinear dynamo models, one can use the kinematic approxima- 
tion to describe the spatial profile of the solution. Certainly, in the nonlinear 
case, spectral parameter 7 differs from eigenvalues (|2"5f) (in the steady non- 
linear regime, 7 has only an imaginary part which determines the period of 
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the solar cycle). However, our analysis of the dynamo wave behavior near 
the pole remains applicable. 
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